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I.*  INTRODUCTION 

The  TROTT1  computer  code  is  being  used  to  simulate  explosive  and  plate 
impact  densifications  of  various  materials. 2  in  the  present  report  the 
effects  of  two  computational  parameters,  cell  size  and  slide  lines,  on  com¬ 
puted  compactions  will  be  discussed. 

TROTT  is  a  Lagrangian  finite  difference  computer  program  for  calculat¬ 
ing  two-dimensional  stress  wave  propagation  through  solid,  porous  or  com¬ 
posite  materials.  The  stress  waves  may  be  caused  by  detonation  of  an 
explosive  or  through  an  impact  in  either  planar  or  axisymmetric  flow.  The 
computations  involve  the  standard  leapfrog  method  of  Von  Neumann  and 
Riohtmyer,3  using  artificial  viscosity  to  smooth  shock  fronts.  For  the 
calculations,  the  material  is  divided  into  small,  quadrilateral  or  tri¬ 
angular,  cells  or  continuum  elements.  The  motion  and  stresses  throughout 
the  material  are  determined  as  a  function  of  time  by  solving  the  mass, 
momentum  and  energy  conservation  relations,  together  with  the  constitutive 
relations  of  the  material. 


II.  COMPUTATIONAL  PROCEDURES 

The  BRL  is  studying  processes  by  which  powders  are  dynamically  com¬ 
pacted  into  solids.  The  experimental  portion  of  these  compaction  studies 
is  supplemented  by  the  use  of  computer  simulations.  Material  variables 
that  are  taken  into  account  in  the  computer  simulations  include  type  and 
amount  of  explosive,  sample  geometry,  container  material,  initial  denr  :ty 
of  powder,  and  temperature.  It  is  essential  that  the  computational  ,  ■  cess 
itself  does  not  influence  the  results  of  the  compaction  simulation.  The 
purpose  of  the  present  report  is  to  assess  the  effects  that  computational 
cell  size  and  slide  lines  have  upon  the  ability  of  the  TROTT  Computer  Code 
to  accurately  simulate  dynamic  compaction  of  various  powders. 

A.  Cell  Size  Effect. 

In  choosing  the  dimensions  of  computational  cells,  too  large  or  too 
small  a  cell  should  be  avoided.  If  too  small  a  cell  is  used  the  computa¬ 
tions  become  excessively  long  and  expensive  whereas  too  coarse  a  cell  size 
can  lead  to  misleading  and  inaccurate  results.  One  way  to  ascertain 
whether  or  not  a  cell  size  is  adequate  is  to  perform  calculations  in  which 


1L.  S  aman  and  D.  R.  Curran,  Contract  Report  ARBRL-CR-00428 ,  "TROTT  Computer 
Program  for  Two -Dimensional  Stress  Wave  Propagation, "  SRI  International , 

333  Ravenswood  Ave.,  Menlo  Park,  CA  94025,  April  1980,  AD-A085766 . 

2G.  L.  Moss,  R.  F.  Benck,  P.  H.  Netherwood,  and  J.  R.  Stratton,  "Dynamic 
Compaction  of  Ceramic  Materials , "  Ballistic  Research  Laboratory ,  Memoran¬ 
dum  Report,  ARBRL-MR-03327 ,  Aberdeen  Proving  Ground,  MD  21005-5066 , 
December  1983.  AD-A137257. 

V.  Von  Neumann  and  R.  D.  Richtmyer ,  "A  Method  for  the  Numerical  Calculation 
of  Hydrodynamic  Shocks,"  J.  Appl.  Phys. ,  Vol.  21,  p.  232  (1950). 
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the  size  of  the  cells  is  varied.  If  the  results  of  the  calculations  are 
essentially  identical  it  can  be  reasonably  assumed  that  the  results  are  not 
dependent  on  cell  dimensions,  at  least  over  the  dimension  interval 
involved.  Cell  size  effects  will  be  determined  for  compactions  accom¬ 
plished  with  impacts  and  explosions. 

1.  Explosive  Compaction. 

Processes  that  involve  detonations  are  simulated  in  TROTT  by  using 
the  EXPLODE  subroutine.  EXPLODE  is  based  upon  Chapman-Jouguet  detonation 
theory  which  incorporates  the  concept  that  transformation  between  chemical 
energy  and  internal  energy  occurs  within  the  shock  front.  (For  further 
discussion  see  any  standard  reference  such  as  Taylor.)* 

Figure  1  shows  a  compaction  configuration  used  in  a  test  of  the 
cell  size  dimensions  on  computed  results.  The  copper  powder  is  60  percent 
of  theoretical  density  (60$T)  and  is  contained  in  a  copper  holder  posi¬ 
tioned  between  two  layers  of  amatol  (0.9*  g/cm^OjtTNT-SO^NH^NO^) .  The  com¬ 
plete  simulation  involved  a  sample  length  of  152.*  mm,  but  only  the  first 
118.8  mm  are  shown  in  Figure  1.  Compaction  is  achieved  by  detonating  the 
amatol  along  the  edge  marked  A-A  in  Figure  1 .  The  simulation  has  been  so 
posed  that  it  results  in  plain  strain  geometry  with  X  and  Y  in  the  plane 
and  Z  in  the  third  (zero  strain)  direction.  Hugoniots  of  solid  and  porous 
copper  as  interpreted  by  Seaman  and  Linde5  were  used  to  characterize  the 
copper.  A  subroutine  POREQST,^  developed  as  a  dynamic  response  model  for 
constitutive  relations  for  porous  materials,  was  used  to  simulate  consoli¬ 
dation  of  the  copper  powder. 

In  Figure  1  the  problem  has  been  divided  into  a  60  (horizontal)  by 
58  (vertical)  grid.  The  effect  of  a  coarser  cell  grid  was  determined  by 
rerunning  the  Figure  1  simulation  using  a  *0  (horizontal)  by  29  (vertical) 
grid.  A  comparison  of  the  two  configurations  is  made  in  Table  1. 


4J.  Taylor ,  "Detonation  in  Condensed  Explosives , "  Clarendon  Press,  Oxford, 
England ,  1952 . 

SL.  Seaman  and  R.  K.  Linde ,  AFWL-TR-68-14Z ,  "Distended  Material  Model 
Development ,  Volume  1,  "  Stanford  Research  Institute ,  Menlo  Dark,  CA 
94025 ,  May  1969 ,  AD-855418. 

6L.  Seaman ,  R.  E,  Tokheim,  and  D.  R.  Curran ,  DNA  Z412F,  " Computational 
Representation  of  Constitutive  Relatione  for  Porous  Material ,  "  Stanford 
Research  Institute ,  ZZZ  Ravenewood  Ave.,  Menlo  Park,  CA  94025,  May  1974, 
AD-A007921. 
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TABLE  I 


Cell  Dimensions  Used  in  Copper  Compaction  Studies 


MATERIAL 

FINE  GRID 

COARSE  GRID 

Cell  dimensions,  mm 

Cell  dimensions,  mm 

Y 

X 

Y 

X 

vertical 

horizontal 

vertical 

horizontal 

amatol 

1.33 

2.58 

2.86 

3.91 

copper 

container 

1 .27 

3.91 

3.18 

3.91 

copper 

powder 

1.27 

2.58 

1 .91 

3.91 

2.  Plate  Impact. 

TROTT  was  used  to  simulate  the  impact  of  a  mild  steel  flyer  plate 
onto  a  mild  steel  cylinder  containing  a  cavity  filled  with  powdered  alumi¬ 
num  oxide,  Al^O^ .  The  velocity  of  the  flyer  plate  at  impact  was  281  m/s. 

A  60  (horizontal)  by  59  (vertical)  rectangular  cell  grid  representation  of 
the  system  at  impact  is  shown  in  Figure  2.  For  this  calculation  an  axi- 
symmetric  geometry  with  X  axial,  Y  radial  and  Z  circumferential  was  used. 
The  axis  of  symmetry  in  Figure  2  is  the  upper  edge  of  the  grid.  To  compute 
the  effect  of  a  coarser  cell  grid  the  compaction  was  rerun  using  cells 
approximately  three  times  as  large  as  those  used  in  Figure  2.  In  both  com¬ 
putations  the  cell  sizes  in  the  flyer  plate,  Al^O^  powder  and  the  portion 

of  the  sample  holder  Detween  flyer  plate  and  A^O^  were  the  3ame  throughout 

each  region.  The  cells  in  the  remainder  of  the  sample  holder  were  kept  at 
constant  width  (Y)  but  the  length  (X)  of  each  succeeding  row  of  cells 
beyond  the  Al^O^  was  increased  by  approximately  12  percent.  This  steady 

increase  in  ce_l  size  allowed  the  existing  experiment  to  be  modeled, 
avoided  awkward  cell  discontinuities  and  limited  computing  time. 

B.  Slide  Line  Effect. 

TROTT  incorporates  a  slide  line  routine  that  allows  pressures  to  be 
transmitted  through  slide  line  boundaries  but  minimizes  distortion  in 
boun^ry  layer  cells.  Minimum  distortions  are  usually  desirable  but  slide 
lines  can  add  an  extra  element  of  artificiality  to  the  results.  One  of  the 
purposes  of  the  current  report  is  to  model  real  phonomena  so  computations 
with  and  without  the  use  of  slide  lines  were  investigated.  Calculations 
were  made  for  an  explosive  compaction  of  powdered  Al^O^  in  a  molybdenum, 

Mo,  container;  specifically  20  mm  of  amatol  was  placed  on  a  6.35  mm  thick 
Mo  plate  placed  on  a  12.5  mm  thick  layer  of  Al^O^  powder.  The  bottom  edge 


FLYER  PLATE 

281  m/s 


Figure  2.  Initial  Conditions  for  Plate  Impact  Experiment. 


of  the  powder  was  constrained  to  act  as  a  rigid  boundary.  The  initial 
computational  cell  representation  for  this  test  is  illustrated  in  Figure  3 
with  the  slide  line  imposed  between  the  amatol  and  Mo  wall.  Computations 
were  made  without  the  slide  line  for  comparison.  A  planar  detonation  was 
initiated  along  the  left  edge  of  the  explosive. 


III.  RESULTS 

The  effects  of  computational  cell  size  and  slide  lines  upon  simula¬ 
tions  of  explosive  and  impact  driven  powder  compactions  were  computed. 

A.  Cell  Size  Effect. 

The  effect  of  increasing  computational  cell  size  upon  the  computed 
density  and/or  pressures  of  powdered  copper  and  Al^O^  were  investigated  for 
both  explosive  and  impact  driven  densifications . 

1.  Explosive  Compaction  of  Copper  Powder. 

The  effects  of  the  detonation  of  amatol  on  the  copper  powder  den¬ 
sity  using  initial  configurations  as  given  in  Figure  1,  were  computed  and 
are  presented  in  Figure  4.  The  density  presented  in  Figure  4  is  that  of  a 
horizontal  slice  of  the  copper  powder  measured  along  the  center  of  the 
copper  powder.  The  area  sampled  is  along  B-B  of  Figure  1.  Figure  4  repre¬ 
sents  the  computed  copper  powder  density,  28.8  ps  after  line  detonation  of 
the  amatol  and  does  not  indicate  the  final  density  of  the  copper.  Calcula¬ 
tions  at  later  times  show  that  the  copper  eventually  attains  a  density  of 
100;ST.  As  Figure  4  illustrates  the  calculated  densities  from  the  two  grids 
are  essentially  equivalent  along  BB. 

2.  Plate  Impact  Compaction  of  AlpO^. 

Figure  5  is  a  plot  of  calculated  Al-,0^  density  as  a  function  of 

time  after  impact  for  a  cell  located  10  mm  from  the  center  line  (axis  of 
symmetry)  and  4  mm  behind  the  rear  side  of  the  cover  plate  (edge  effects  at 
this  location  should  be  minimal).  Figure  5  includes  data  generated  from 
both  fine  and  coarse  cell  grids.  Figure  5  indicates  that  there  appears  to 
be  little  difference  between  densities  generated  through  the  use  of  the  two 
grid  sizes. 

B.  Slide  Line  Effect. 

Figure  6  shows  the  computed  geometrical  representation  of  a  sweeping 
explosive  lensification  of  AlpO^,  42.4  ps  after  detonation  using  the  slide 

line  routine.  A  similar  curve  at  41.9  ps  for  a  calculation  without  using 
slide  lines  is  presented  in  Figure  7.  Computed  copper  densities  with  and 
without  slide  lines  are  presented  in  Figure  8.  These  calculations  are  at 
33. ’0  ns  after  detonation  initiation  and  are  an  average  of  the  cell  densi¬ 
ties  at  each  vertical  location.  As  is  illustrated  in  Figure  8  the  powder 
is  fully  densified  in  the  first  two  cm  and  doesn't  appear  to  be  overly  sen¬ 
sitive  to  the  use  of  slide  lines  even  though  the  expansion  of  the  explosive 
products  is  noticeably  different. 
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IV.  CONCLUSIONS 


TROTT  computer  simulations  indicate  that  under  the  experimental  condi 
•,i:ns  discussed  herein  for  the  explosive  compaction  of  powdered  copper  and 
plate  impact  densif ications  of  Al^O  a  30  by  MO  grid  of  quadrilateral 
ceils  is  adequate.  J 

TROTT  calculations  also  indicate  that  in  the  simulation  of  a  sweeping 
dunsification  of  AljO^  the  presence  or  absence  of  a  slide  line  does  not 

aocreciably  affect  the  final  calculated  A1.0.  density. 


WITHOUT  SLIDE  LINE 


Figure  8.  Density  of  Copper  Powder  as  a  Function  of  Distance  from  Solid 
Copper  Interface,  33.1  us  after  Detonation  Initiation. 


Figure  7.  Amatol  -Mo -A  l^,  41.9  ys  after  Detonat 
Initiation  -  Without  Slide  Line. 
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Figure  5.  Density  of  Alo03  Powdfc)"  as  a  Function  of  Time  after 
Plate  Impact  for  Both  Fine  and  Coarse  Grids. 
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